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Distributions of charge near charged impurities in Weyl semimetals are considered with the help 
of relativistic Thomas- Fermi method in full analogy with the solutions previously found in QED. 
Screening and anti-screening, zero charge and asymptotic freedom solutions appearing in different 
physical situations are discussed. 
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Introduction.— Existence of the Weyl semimetals, i.e. 
the materials with the points in Brillouin zone, where the 
completely filled valence and completely empty conduc- 
tion bands meet with linear dispersion law, u> = v-pp, has 
been predicted ini. The order of value estimate of the 
Fermi velocity is vy ~ 10~ 2 c, where c is the speed of light 
in vacuum. Systems with relativistic dispersion law are 
likely to be realized in some doped silver chalcogenides^, 
pyrochlore iridates^, and in topological insulator mul- 
tilayer structures^. Weyl semimetals are 3-dimcnsional 
analogs of recently discovered graphene^, where the en- 
ergy of excitations is also approximately presented by 
the linear function of the momentum but the electron 
subsystem is two dimensional one whereas the photon 
subsystem remains three dimensional. Even though the 
mass of excitations m = for ideal graphene and Weyl 
semimetals without interactions, m ^ can be induced 
in many ways resulting in the gapped dispersion relation^ 



= p 2 Vp + m 2 Vp. To be specific and in order to avoid 



discussion of infrared divergence we consider below the 
gapped case. In difference with a weak value of the fine 
structure constant in QED, cvqed = e 2 /hc = 1/137, 
the coupling in Weyl semimetals and in graphene is 
rather strong, a = e 2 /hvF ~ 1. The effective coupling 
a e { = e 2 /hv-peo, where eo is the dielectric permitivity of 
the substance, can be as J 1 as < 1 depending on the 
substance and both weak and strong coupling regimes 
are experimentally accessible. Thus Weyl semimetals and 
infinite stack of graphene layers make it possible to ex- 
perimentally study various effects have been considered 
in 3+1 quantum electrodynamics (QED) for weak and 
effectively strong couplings. 

Here we study the screening problem in Weyl semimet- 
als, infinite stack of graphene layers and 3+1 QED. For 
that we use the relativistic Thomas-Fermi approach de- 
veloped in QED and applied to find charge distribu- 
tions near the extended charge source 7 - and the point- 
like charge source^. Recently RefJ^ argued for the zero 
charge effect for the impurity screening in Weyl semimet- 
als in the limit of vanishing impurity size. We argue for 
another solutions. 

Single-particle Coulomb problem. Falling to 
the center.— The energy of the electron levels in the 
Coulomb field V = —Zoe 2 /r of a point charge Zq is found 
from the solution of the Dirac equation, seeii, 



Here K = + (j + 1/2), j is the total angular momen- 
tum, n r is the radial quantum number, n r = 0,1,2... 
for k < and n r = 1,2,... for k > 0, Qq — ZoOqed 
in QED case. Eq. ([T]) is symmetric under the simulta- 
neous replacements ui — > —u, Zq — > — Zq. Thus, even 
for small Zq > moreover the electron levels follow- 
ing from the upper continuum with w' 1 ' < 



' mc 2 there 



are levels with 



,(2) 



A 1 ) 



> —mc 2 following from 



the lower continuum. According to the traditional inter- 
pretation supported by experiments electron states with 



,(2) 



> 



are interpreted as positron states for Zq < 



with U> + = -wf = We . 

The energy of the ground state electron level (n r = 0, 
k = —1) appearing from the upper continuum decreases 
with increase of Zq > 0, luq = mc 2 y/l — Qq, and vanishes 
for Z c = 1/cvqed- For Zq > I/chqed the single particle 
problem for a point nucleus loses its meaning because 
of the collapse to the center. Eq. (JTJ) is also valid for 
Weyl semimetals after replacement c — > v-p. Thus for 
semimetals the falling to the center in the single particle 
problem may occur already for Zq > Z c = 1/a ~ 1 . For 
masslcss particles in the Coulomb field the bound state 
appears for Zq > Z c and the critical charge for the falling 
to the center is the same as for the massive particles. 

Finite size nucleus.— To resolve the problem Ref.— 
suggested to consider the central charge as smeared out 
over a sphere of the finite radius R. For the nucleus of 
the finite radius the same as that for ordinary atomic 
nuclei, i.e. tqA 1 ^ (for tq ~ 1.2 fm and atomic number 
A ~ IZq), the electron ground state level in the QED 
case crosses the upper boundary of the lower continuum 
at nuclear charge Zq ~ 170, see^. The energy of the level 
with id < —mc 2 (Z > Z c and/or R < R c ) acquires a neg- 
ative imaginary parfci^ cj = Rew + ilmuj, Imw < 0, corre- 
sponding to a decay of quasi-stationary state. Two elec- 
trons penetrate from the lower continuum to the ground 
state level in upper continuum that corresponds to tun- 
neling of two positrons (the holes for Weyl semimetals) 
to infinity. For fixed Zq > 1/a the energy of the ground 
state electron level continues to drop with decrease of 
R and crosses the upper boundary of the lower con- 



tinuum (valence band, w_ 



for Weyl semimet- 



\k 2 



fo) 1/2 



(1) 



als). It occurs at the impurity charge Zq — > 1/a + 
(l/a c f if we included screening) for R — > 0, se o 15 ' 16 . As 
was believed during long time, at least for Zq close to 



2 



Z c , owing to Imw < for the overcritical level and the 
Pauli principle, smearing of the charge source even for 
R — >• allows to solve the problem. However in real- 
ity it is not the case, since at small distances and in 
a strong field there arises a strong polarization of the 
vacuum and the bare Coulomb potential V = —Zoe 2 /r 
should be replaced by V — — ZqC 2 j t{y)r , where e(r) is 
the dielectric pcrmitivity. The later quantity decreases 
with decrease of r that effectively corresponds to increase 
of the charge at small distances. On the other hand, 
with increase of Zq and decrease of R many levels cross 
the boundary of the lower continuum (valence band). 
The level density, p(cu < — mVp) = dn/duj ~ \Jy for 
< y = R\/oj 2 — mPv^/y/Q 2 , — 1 < 1, is crowding to- 
ward the boundary of the valence bancU^. After tunnel- 
ing electrons occupy all available levels with ui < —mc 2 
(vacuum shell) in QED and u> < —mvp in case of Weyl 
semimetals. Thus the problems of the screening of the 
point-like charge with Zq > 1 /a, as well as of an extended 
charge source with Zq 1/at, are actually many-particle 
problems. 

Many-particle problem. Electron condensation 
in upper continuum.— Consider a spherically symmet- 
ric charged impurity with the number density rii mp (r) 
embedded into a semimetal. Let V = — eAo is the self- 
consistent potential for the electron. Since many elec- 
tron states condense (in the sense that electrons in accor- 
dance with the Pauli principle occupy all available energy 
levels) in presence of impurity with a sufficiently large 
charge, they have in average large angular momenta. 
Thereby spin effects (except Pauli principle) are inessen- 
tial and the electron (Fermi) momentum can be found 
from the Klein- Gordon equation 



h 2 Ail> +p 2 {r)i> = 0, 



p 2 (r) = {LO-V) 2 /v 2 F -m 2 vl.{2) 



The quasiparticlc approximation holds at <C 1 

that for the Coulomb field, V = —Z e 2 /r, yields criterion 
Zqol 3> 1. However note that quasiclassical approxima- 
tion reproduces exact Eq. (JXJ) for purely Coulomb field for 
any Zq and has a good numerical accuracy also for finite 
size nucleus with Zq > l/a G f, see^. E.g., with percent- 
age accuracy it allows to get expression for single-particle 
energy of electron levels below the upper boundary of the 
valence band, 



-mv\ 



R c /R-( K -l/2)/2Q 2 
1 - (« - l/2)/2Q§ ' 



(3) 



with Qq = Zoa e f, and critical radius R c = R c (n r , n, Zq). 
The classically accessible regions p 2 (r) > correspond 
to the upper and lower continua (conduction and valence 
bands), being curved in presence of the field and the re- 
gion p 2 {r) < corresponds to the gap. The electron lev- 
els are in the upper continuum, if they drive to to = mc 2 
(uj = mvp) with switching of the potential V < 0, and 
they would be in the lower continuum, if they reached 
oj = —mc 2 (oj — —mvp) with switching of the poten- 
tial. As we mentioned, in the later case following the 



traditional interpretation after the replacement lj — s- — uj, 
k — y — k the levels are interpreted not as electron levels in 
the lower continuum but as levels of positrons/holes (in 
the repulsive potential — V for the electron) in the upper 
continuum. We return to this interpretation below. 

Within the relativistic Thomas-Fermi approximation 
the chemical potential u;_ is determined by the last level 
filled by the electrons —mv\ < cj_ < mv\. If there are 
no free electrons in the system, in the field V only the 
levels of the vacuum shell u < w_ = —mv 2 , are occupied 
by electrons penetrated from the valence band, that cor- 
responds to p 2 _(r) = (cj_ — V) 2 /vp — m 2 Vp. The holes 
(positrons in QED) go off to infinity. If there is sufficient 
amount of electrons in the conduction band (case of the 
neutral atom in QED), then u;_ = mv 2 ,. 

The number density of filled electron states in upper 
continuum is given bySS 



n e (r) = gp 3 _(6ir 2 h 3 ) 



3\-l 



(4) 



where g is the degeneracy factor, g = 24 for py- 
rochlorc iridates^ and g — 2 in a topological insula- 
tor multilayer— and for electrons in QED. With the 
Coulomb law the number of electrons in electron con- 
densate N e = J R 4:Trr 2 drn e oc — lni? diverges at R — > 
showing that "electron condensate" should significantly 
modify the Coulomb law at small distances. 

The relativistic Thomas-Fermi equation for the de- 
scription of the charge distribution is as follows, cf.— : 



V(e(r)W(r)) = 4^e 2 [n imp (r) - n e (r)}, 



with the boundary conditions at the ion radius r^. 

V(ri) = uj- — mv 2 , > — 2mvp , 
V'(n) = (Z(r)e 2 /e(r)ry r ^ + . 



(5) 



(6) 



The first condition follows from the requirement n e (ri) 
0. In case of the empty conduction band cj_ = 



F • 



and for the metal u>- = mi) 2 ?, r% = oo. At distances, 
where the screening is most effective \ V\ 3> mv 2 ? and Eq. 
([5]) simplifies as 

V(eW) =4^e 2 [n irap + 5 ^ 3 /(67r 2 ^)], V < 0. (7) 

1. Limit of strong coupling. Let the bare charge 
distribution in impurity is rii mp = uqO(R — r), uq — 
Z /(4:TtR 3 /3) = const, 9(x) = I for x > and for 
x < is the step function. Assume R^$> I. Then at typ- 
ical distances x = (i — R)/l ~ 1 the geometry is reduced 
to the flat one. Also assume that the dielectric permi- 
tivity e varies at distances l e 3> I, where I is the typical 
size of the change of the electric field. Then at distances 
~ I of our interest in the given case one may not vary e. 
Replacing V = — Vqx(x), in Eq. we find, cf.— : 

X = x 3 -0(-x), V = hv F (6ir 2 n /g) 1 / 3 , (8) 



e a e f V 32 ^"^ 



1/3 



eR 2 



3 2 /3 eo(A Z 2 )i/3 ' 
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where A = a£ f . Inequality Vb ^> mv\ yields n ^> 
g(mvF/h) 3 /(6iT 2 ). The used condition R^> I is satisfied 
for (AZ^) 1 / 6 > (e/e ) 1/2 for arbitrary i?, even for R 0. 



Solutions of (TST) are as follows 



X(x) = 1 - 



l + 2- 1 /2 s h(a-a;/y3) 
X(x) = v^/(x + 6), a; > 0. 



, x < , (9) 
(10) 



Matching of these solutions with derivatives at x = 
yields values of constants: sha = ll\/2 and b = Ay/2/3. 
Note that the solution (flU)) would have a pole at r = r polc , 
if it were analytically continued in the region x < 0. In 
reality the pole does not manifest since R > r po \ c . The 
solutions ©, (Unj) demonstrate that for (AZ 2 ) 1 / 6 > 1 
(at e(i?) ~ eo)j i- e - i n the strong coupling limit the 
impurity interior is neutral, whereas the charge is re- 
pelled to a narrow (of the length / <C R) layer near the 
droplet surface. The total charge Z{R) situated inside 
the impurity is found with the help of the Gauss the- 
orem, ciX, Z(R)/Z = 3 8 / 3 ■ 2~ 9 / 2 (AZ2)" 1 /6. Solution 
(|TU)l becomes invalid at r ^> R when the charge Z(r) 
decreases up to values Z(r) ~ 1/A 1 / 2 . For still larger r 
the system is described in the effectively weak coupling 
regime: AZ 2 (r) <C 1- In case w_ = mv 2 , the charge tends 
to zero but since this residual screening occurs at r 7S> R, 
the solution does not manifest the Landau zero charge 
effect, i.e. full screening of any bare charge at r ~ R, for 
R —> 0. For the gapless case (m = 0) the screening stops 
for Z(r) < l/a c {, since for such Z there are already no 
electron levels in the Coulomb field. 

2. Limit of weak coupling. Existence of a minimal 
radius. Introduce convenient variables 

V = Qi(r)/r = -Z{r)e 2 /r, t = -\n(r/a). (11) 



Although a is an arbitrary constant it is convenient to 
make a specific choice, e.g. a ~ r^ 7 -^, a = fi / (mvp)^, or 

a = m. 

In the new variables Eq. ([5]) becomes for r > R 



d c{t)Z(t) 
dt 



C-0 



d_ 

dt 



e(<) dZ(t) 



= \z 3 (t)c*e(0(i2) 



i 



2co^eoae 
Z(t)e 2 



Eq dt 

1 (w 2 — m 2 Vp)e 2 a 2 e 
" + Z 2 (t)e 4 



2.1. Approximation e = eo = const. The spatial de- 
pendence of e can be disregarded for rather extended 
charged objects. With e = e , Eq. (fT2j) reduces to 

+ = AZ3(t)^6(0, r>R. (13) 

In QED the solution of this equation has been found in 7,8 
and matched with the Coulomb law for r = r^. 

As has been shown in 8 -, Eq. (fT3|) has the pole solution 



for tp i c - t < 1. In the strong coupling limit the first 
term transforms to the solution (|10p after change of vari- 
ables. Near the pole the second term in the l.h.s. of Eq. 
(fT3f is the dominant one, i.e. |^|| <C whereas the 

first term determines the correction in brackets of (|14|l . 
The value t po ie( A) can be obtained by numerical solu- 
tion of Eq. (fT2"j) with the boundary condition ©. For 
oj_ = -rai'p, v-p = c, 6q = 1, a = 7-j this solution was 
found in 8 -. In the weak coupling limit, (AZ(r.i) 2 ) 1 / 6 <C 1, 
the value t po \ c = 1/8/z — InD, where lnD(/i) is a smooth 
function of fi = \Z 2 (t = 0)/(2g), (e.g., D ~ 0.2/i 1 / 2 for 
0.4 < /i < 1), see Fig. 4 in 8 -). This result approximately 
(due to a smooth logarithmic dependence of t po \ on D) 



holds also for 



F > 



«f 7^ c. Thus in the weak 



coupling limit, for any value of the charge at large dis- 
tances Z(r > n) = Zoo the solution Eq. (fT3|) matched 
with the purely Coulomb law for r > exists only if the 
impurity has the radius R > r po \ c . 

In the limit > |0|, for r < r<, i.e. for f ~ 1, 
Eq. (fT3")) also has analytic solution: 



Z(i) = Z(i = 0)/i/l - 2XZ 2 (t = 0)t, 
where Z(t = 0) = Z(r = a). Comparing values 



(15) 



\ — \ 



and 



dT 1 



we see that solution (fT5|) holds only for \Z 2 (t 



0)/(l - 2AZ 2 (t = 0)t) < 1. As we see, for a = n (then 
t > 0) this inequality holds only in the weak coupling 
limit XZ 2 (t = 0) < 1 and at distances when l-2AZ 2 (i = 
0)< > AZ 2 (i = 0). Taking a = n, r = R < n in dUD 
we recover the relation between the observable charge 
= Z(r > Ti) and the charge Z(r = R) at the surface 
of impurity 

Zoo ^ Z(i?)[l - AZ 2 (i?)ln^], at AZ 2 (i?)ln^ < 1(16) 

it it 



Also on example of Eq. (|T3| Rcf.— demonstrated prob- 
lems with the usage of the leading logarithmic approx- 
imation (LLA), being often exploited beyond its region 
of validity in problems of quantum field theory. Expend 
Z{t) in the series of the parameter /Lt, Z(t) = ^2 n Z n /j, n , 
with Z„ = C n t n + Oit 71 ' 1 ), for t -> oo. For the leading 
coefficients C n we get the recurrence relations 



Co — 2 , C„ — — 



(17) 



which yield C„ = 2" +1 (2?i)!/(n!) 2 . Summation in LLA 
recovers solution (fT5|l . Note that (fT5|) has the square-root 



singularity for £ — > <lla — iTz^Xt^o)' wncreas the real 
solution has pole for t — > t po i < fLLA- This is because 
the LLA becomes invalid in the region where 1 — XZ 2 (t = 
0)t < XZ 2 (t — 0). In this region sub- leading corrections 
^n-j-n-i om j^ted in the derivation should be incorporated. 

On the other hand, solution (fT5|) coincides with the 
one obtained ir>iS for a — R (t < 0): 



Z(t) 



^pole 



1+ g(tpole-f) + 



(14) 



Z(r> J R) = Z(i?)A/l 



2XZ 2 (R) In-. 



(18) 
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Now the criterion |^|| > 

2\Z 2 (R)\nj i ) < 1, i.e. for all t < 0, but only for 
\Z 2 (R) <C 1, i.e. in the weak coupling limit. Con- 
trary, RefJ^ used solution (fT8|) for any r > R at arbi- 
trary small R in the strong coupling limit XZ 2 (R) ^> 1 
pointing on the zero charge behavior of the charge dis- 
tribution, i.e. dropping of the charge Z(r) to zero in the 
narrow vicinity of r = R. However (|15[) does not hold 
in this strong coupling limit at r ~ R. Also, for gapless 
substances, as considered ini£, the relativistic Thomas- 
Fermi approximation holds only for Z(r)a e { 3> 1 and 
thus solution (fT8|) is not correct for r S> R when Z(r) 
becomes < l/a c f- Thus we do not support conclusion of 
Ref.— on the zero charge behavior of the solution of the 
relativistic Thomas-Fermi equation for the point charge 
source. 

Concluding above discussion, the receipt to consider 
the point source of charge as smeared one for R —> to 
resolve the problem of the falling to the center, which 
worked in the single-particle problem, does not work in 
the many-particle problem in the weak coupling limit. 
Repeating the result^ we stress that for e = const and 
XZq <C 1 relevant solutions of Eq. (fT2"j) reproduce the 
Coulomb law at large distances, but only for extended 
charged droplets with R > r po i e - 

2.2. Taking into account spatial dispersion of e. The 
dielectric permitivity of the 3+1 QED vacuum in elec- 
tric field can be found with the help of the interpolation 
formula, valid with logarithmic accuracj 



l dp 
for 



holds for XZ 2 (R)/(1 + describing both large and small distances is as follows 



17. 



/ \ , «QED , 

eqed \T) = 1 q — In 

07T 



eETi 



(19) 



where eE — —V'(r) = Q{r)/r 2 is the electric field ten- 
sion. From this equation in case of the strong homoge- 
neous electric field we recover the result of Heisenberg 
and Euler—. In case of the strong inhomogeneous elec- 
tric field, for Q(r) > 1, Eq. ([19]) yields the solution^ 
that generalizes the Heisenberg and Euler result. In case 
of a weak field we recover the Uehling and Scrber cor- 
rection at r < h/mc (see second term (fT9)) '). Also Eq. 
(|T9| correctly transforms to the result derived in LLAii. 
Note that the latter result is derived at the condition 
ln(fi,/rmc) 3> 1 being formally valid as for e > as for 
e < (for |e| > a QE D/37r). 

At small distances the field becomes strong indepen- 
dently of whether the observed charge Z^ is large or 
small. Indeed, in absence of the electron condensation 
Z(r) = Z oc /e(r) and dielectric permitivity decreases as 
r decreases. Thus, at small distances any case we have 
Z(r) > 1/ctQED and the electron condensation occurs at 
levels in the upper continuum, which have crossed the 
boundary w_ = —mc 2 . 

In matter at large distances r>I e one has e ~ const = 
e . But at much smaller distances and/or in presence of 
very strong electric field the screening should occur sim- 
ilar to that in vacuum. Thus the interpolation formula 



<r) = eo 1 



ga c f 
Gtt 



In 



eEh 



(20) 



Assume 



~3t~\ ^ lei ( e w) I- Tnen disregarding sec- 
ond term in the l.h.s. of Eq. (fT2|. setting £ = 1 in the 
r.h.s. and using that at the above condition the charge 
is a smooth function, Q ~ Qi, we arrive at the solution 



Z 2 (r) = (-2a 2 cf + Ce 2 y 



(21) 



Constant C can be found from interpolation of the so- 
lution to the Coulomb law at large distances, i.e. from 
that e — > eo, Z — > Z^ for r — > oo. Thus we obtain 



Z 2 (r) 



(l + 2a 2 cl Z 2 oo )e 2 (r)/e 2 



(22) 



This equation has the inflection point at t = t- m & at which 
dZ/dt = oo. However already before reaching this point, 
i.e. for e - v /2a^ f Z^ 3 /(l + 2a 2 { Z 2 <3 ) - ga ci /6n, inequal- 
ity 



I d{eZ) 



i ^> \4i (e^) | becomes invalid. Z(r) contin- 
ues to grow with decrease of r. For smaller r inequality 

I djeZ) I „ / dZ 
dt I ^ I dt V fc dt 



<C 1^7 (e^T-) | is fulfilled and the pole- like solution 




is generated. Using that e is still a smooth function of t, 
except very near the pole, we obtain 



(23) 



similar to (|T4j) . This solution ceases to be useful only 
for e < ga c f/6n, since the condition of smoothness of 
variation of e becomes invalid. In this region, before oc- 
curring of the pole but very near the pole, the solution 
(|23| has a point of inflection t = tj n fl. Thus Eq. (fl~2]) has 
a solution that falls off with increase of r and reduces to 
the Coulomb solution only if the radius of the source is 
R > r infi > r poic- This is how the falling to the center 
manifests in the many-particle problem provided spatial 
dispersion of e and electron condensation on the levels 
in the conduction band (upper continuum) are included. 
We see that the initial problem of the falling to the centre 
became more severe than in the single particle case where 
any smearing of the charge source (even for R — > f)2i) was 
sufficient to overcome difficulties. 

Impurity of an arbitrary small size. Hypothe- 
sis of electron condensation in lower continuum. 
In Ref.— a solution of the problem for the charge distri- 
bution near a source of the radius R < r po \ c was pro- 
posed in case of QED. One observes that the dielectric 
permitivity given by Eqs. (fT9|) and (|20| becomes nega- 
tive for r < r e , r e is the point where e = 0. Thereby 
attraction in the original potential (for Zq > 0) is re- 
placed by an effective repulsion to the electron at small 
distances, V > 0. On the other hand, in the repulsive 
potential cither positrons (holes) are accumulated in the 
upper continuum (this case after replacement Z — > — Z 
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yields the same solution, as we have found for electrons 
in the attractive field, being valid only for r > r po i c ) or 
the electrons are condensed right in the lower continuum 
(without any tunneling) as it is allowed by the symmetry 
uj — > — w, V — > — V of Eq. (Q}. Thereby now consider 
possibility of the electron condensation on the levels in 
the lower continuum in the repulsive potential. 
The charge distribution is then determined by 

V(eW) = 47re 2 [7i imp - gV 3 / (6Tr 2 h 3 v^)}, V > 0. (24) 

For a point charge Zq we have ni mp = Zq5(t). Notice 
change of the sign in the second term in the r.h.s. com- 
pared to that in Eq. ([7]) since this term corresponds now 
to electrons condensed on levels of the lower continuum 
in the field V > 0. For r > Eq. reduces to 



at eo 



Assuming | ' | 3> |^ (e^jf) | we now obtain 
Z 2 (r) = (2a c 2 f + C l£ 2 )- 1 . 



(26) 



Constant C\ can be found from the condition Z 2 (r 
0) — > Zq since then e 2 — > oo. Thus we obtain 



Z\r) 



7 2 



2a 2 f Z 2 



(27) 



The solution is similar to that found in QCD with taking 
into account of the quark condensation in gluo-electric 
fielcU^. A similarity and difference between QED and 
QCD solutions were analyzed in RefJ^. Solution ([27]) 
corresponds to asymptotically free regime at extremely 



small distances Z(r -> R . 0) -> -Z /\e(R)\ ~> 0. No- 
tice difference in the sign of Z(r) and Zq. Z 2 (r) grows 
with r and reaches maximum at e = (for r = r e ). 
Then Z 2 {r) decreases with subsequent increase of r. At 
r > r e , e -> e and Z{r) -> Z M = -Z / V 2a of Z o + e o- 
For Zo <C 1/oief, ^oo — — ^o/ £ o- In QED such a so- 
lution proves to be consistent with the renormalization 
relation between the bare and physical charges that ar- 
gues in favor of consistency of QED as theory with point 
interaction. It would be very interesting to experimen- 
tally check this peculiar possibility, e.g. measuring the 
field of a nucleus embedded in the Weyl semimetal. 

3/2 

Concluding, in the strong coupling limit, Zq ^> a el , 
owing to the screening the interior of the impurity em- 
bedded in a Weyl semimetal becomes electrically neutral 
even for R — > and the charge is repelled in a narrow 
layer near the surface, Ar <IC R. With taking into ac- 
count of the electron condensation on the levels in the 
upper continuum (conduction band) and e(r) > we also 
found the charge distribution for impurity with a radius 
R > r po \ c in the weak coupling limit, a e{ <C Zq <C a f . 
For the charge source with the radius R < r po \ e there also 
exists a solution demonstrating asymptotic freedom, pro- 
vided negativeness of the dielectric permitivity at small 
distances and the electron condensation on the levels 
in the lover continuum (in valence band). Thus Weyl 
scmimetals give intriguing possibility to experimentally 
check validity of these solutions, in particular to verify 
the hypothesis of the electron condensation on the levels 
of the lower continuum, possibilities of e < and asymp- 
totic freedom at small distances. 
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